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Abstract. Phase transitions in polar binary mixtures can be drastically altered 
even by a small amount of salt. This is because the preferential solvation strongly 
depends on the ambient composition. Together with a summary of our research in this 
problem, we present some detailed results on the role of antagonistic salt composed 
of hydrophilic and hydrophobic ions. These ions tend to segregate at liquid-liquid 
interfaces and selectively couple to water-rich and oil-rich composition fluctuations, 
leading to mcsophasc formation. In our two-dimensional simulation, the corascning of 
the domain structures can be stopped or slowed down, depending on the interaction 
parameter (or the temperature) and the salt density. We realize stripe patterns at the 
critical composition and droplet patterns at off-critical compositions. In the latter case, 
charged droplets emerge with considerable size dispersity in a percolated region. We 
also give the structure factors among the ions, accounting for the Coulomb interaction 
and the solvation interaction mediated by the composition fluctuations. 
(Some figures in this article are in colour only in the electronic version) 
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1. Introduction 



In many soft materials, phase transitions and structure formations occur in the presence 
of charged objects including small ions, charged colloids, and polyelectrolytes in polar 
fluids (mostly in water) [l]-[3]. Much attention has been paid to the consequences of 
the long-range Coulomb interaction, but not enough attention has yet been paid to 
the solvation interaction among charged objects and polar molecules |4]-[7]. In the 
theoretical literature in physics, the solvation interaction has been mostly neglected 
in analytic theories and numerical simulations. In real aqueous systems, a number of 
molecular interactions are operative at short distances under the influence of the long- 
range electric potential |4]. In particular, the hydrogen bonding is of primary importance 
in the phase behavior in aqueous systems, as studied in neutral polymer solutions js]. 
However, small ions like Na"*" form a solvation shell composed of several water molecules 
and are hydrophilic, around which the energy of the ion-dipole interaction typically 
exceeds the energy of the surrounding hydrogen bonds. On the contrary, some neutral or 
ionized particles dislike to be in contact with water molecules and become hydrophobic. 
Hydrophilic and hydrophobic particles deform the surrounding hydrogen bonding in 
different manners |4]. Thus the solvation is highly complex in aqueous mixtures. In this 
paper, we will show that the preferential (or selective) solvation of ions in a mixture 
solvent can drastically influence the phase transition behavior. 

Nabutovskii et al. ^ found a possibility of mesophases in electrolytes in the presence 
of a coupling between the composition and the charge density in the free energy. In our 



theory of electrolyte 10,11 , such a coupling originates from the solvation in polar fluids 



and its magnitude is exceedingly high in many real systems. Recently, including the 
solvation effect in mixture solvents, several theoretical groups have examined the ion 
effects in electrolytes 10 19 , polyelctrolytes 20 ,21 , and ionic surfactants 22 using 



phenomenological models. A review on the static properties in this line was presented 



in Ref.[23]. In the dynamics, a number of problems still remain not well studied 16,19 



On the other hand, a large number of microscopic simulations have been performed. 



for example, to calculate the ion distributions near water-air interfaces 24 and solid 



surfaces 25,26 



In this paper, we are particularly interested in statics and dynamics of phase 
transitions in aqueous mixtures with antagonistic salt, which is composed of hydrophilic 
and hydrophobic ions. While the surface tension of liquid-liquid interfaces is known to 



be slightly increased by hydrophilic ion pairs 27,28 , it can be dramatically decreased 



by antagonistic salt as observed 29 , 30 . Hydrophilic and hydrophobic ions segregate 
at liquid-liquid interfaces on the scale of the Debye screening length k"^, resulting in a 
large electric double layer [TT|[l2|[l6] . They also interact differently with water-rich and 
oil-rich composition fluctuations. Mesophases (charge density waves) then appear near 
the solvent criticality for sufficiently large solvation asymmetry between the cations and 
the anions. 

We mention recent experiments where antagonistic salt was used, (i) Luo et al 
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observed a unique interfacial ion distribution by x-ray reflectivity experiment ISOl. (ii) 



Sadakae et al. 31 added a small amount of sodium tetraphenylborate NaBPli4 to a 
near-critical mixture of D2O and tri-methylpyridine (3MP). This salt dissociates into 
liydrophilic Na"*" and hydrophobic BPhJ. The latter anion consists of four phenyl rings 
bonded to an ionized boron. They found a peak at an intermediate wave number 
Q'm(~ 0.1 ~ k) in the intensity of small-angle neutron scattering. The peak 

height was much enhanced with formation of periodic structures. Moreover, they 
observed multi-lamellar (onion) structures at small volume fractions of 3MP (in D2O- 
rich solvent) far from the critical point |32|, where BPhJ and solvating 3MP form 
charged lamellae. These findings demonstrate very strong hydrophobicity of BPh^ , 
which was effective close to the solvent criticality j3l] and far from it p2l. (iii) Another 
intersting phenomenon is spontaneous emulsification (formation of small water droplets) 



at a water-nitrobenzene(NB) interface 33,34. It was observed when a large water 



droplet was pushed into a cell containing NB and antagonistic salt (tetraalkylammonium 
chloride). This instability was caused by ion transport through the interface. 

The organization of this paper is as follows. In Section 2, we will present a summary 
of the solvation effects, particulaly focusing on the preferential solvation. In Section 3, 
the Ginzburg-Landau model and dynamic equations for such systems will be given. As a 
new result, the surface tension of binary mixtures containing an antagonistic salt will be 
shown to decrease to zero with increasing its density. In Section 4, we will numerically 
examine the mesophase formation with an antagonistic salt in two dimensions (2D). 

2. Solvation effects 

This section is an introduction of the preferential solvation. It provides a starting point 
of our theory in Section 3. 

2.1. Solvation shell and Born theory 

In polar liquids, a number of poalr molecules are attached to a hydrophilic ion such 
as Na"*", Li"*", or Ca^"*" to form a solvation shell due to the ion-dipole interaction [4,35 



(see figure 1). On the other hand, hydrophobic particles interact with water molecules 
repulsively, so they tend to form aggregates in water and are more soluble in oil than in 
water. We may consider mixture solvents composed of polar and less polar components, 
where hydrophilic (hydrophobic) ions are surrounded by molecules of the more (less) 
polar component [5|[7] (see the right panel of figure 1). That is, solvation occurs 
preferentially or selectively in mixture solvents. There can be a variety of mixture 
solvents and, in this paper, we call the two components simply as water and oil. Note 
that solvation is also called hydration for water (or heavy water) and for aqueous 
mixtures containing water. 

In this paper, we introduce a solvation chemical potential /isoi(0) each 
hydrophilic or hydrophobic charged particle, where i represents the particle species and 
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(a) hydration in water (b) hydration in water-nitrobenzene 



Figure 1. Illustration of hydration of Na"*" surrounded by a shell composed of water 
molecules in (a) pure water and (b) water-nitrobenzene. 



(j) is the water composition. It is a thermodynamic quantity obtained after the statistical 
average over the thermal fluctuations. With formation of a well-defined solvation shell, 
the typical magnitude of /ig^j much exceeds the thermal energy ksT . In his original 
theory, Born took into account the polarization of polar molecules around a spherical ion 
using a continuum theory of electrostatics 36 . It is the space integral of the electrostatic 
energy ^^^^/Stt, where E = — V$ in terms of the potential $. Here $ = Zic/er depends 
on the distance r from the ion center with Z^e being the charge. The classic Born formula 
thus obtained is written as 

(/i^„i)Born = Zy/2Bl^e = kBTZ^iB/2Rlr„ (2.1) 

which is applicable for hydrophilic small ions. Here the linear dielectric constant e is 
used, though the electric field in the vicinity of a small ion is strong beyond the linear 
response level (leading to dielectric saturation). The space integral is in the region 
-^ion < ^ < oo, where the lower cut-off Rl^^ is called the Born radius of order lA for 
small metallic ions |4 ,37 . In terms of the Bjerrum length £b = e'^/skBT, the Born form 
(Aisoi)Born itsclf Considerably exceeds the thermal energy fc^T for < i^. Here is 7A 
for room-temperature water and is longer for less polar solvents. The Born formula is in 
accord with the fact that the solvation is stronger for smaller ions. It also indicates that 
the solvation is much intensified for multivalent ions such as Ca^"*" or For mixture 

solvents, the dielectric constant e = e{(f)) depends on the composition (p of the more 
polar component A |38|, changing from of the less polar component B to of the 
more polar component A. Since the Born expression is inversely proportional to e{(j)), 
it indicates strong preferential solvation. However, the Born formula is known to be 
very crude neglecting the formation of the shell structure, the density and composition 
changes (electrostriction), and the nonlinear dielectric effect. Moreover, the solvation 
should be strongly coupled with the surrounding hydrogen bonds in aqueous mixtures. 

We give a simple argument on the condition of the formation of a well-defined 
solvation shell composed of several polar molecules. Let the solvent molecular size 

1/3 . o 

a = Vq be longer than the ion radius ~ Rl^^ (where a ~ 3A for water). For single 
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component polar liquids, we compare the electrostatic energy in the region r a and 
the thermal energy /c^T to find the criterion, 

Ao = Z'^e^/STrekBTa > 1. (2.2) 



For mixture solvents, a phenomenological theory yields the criterion 10 , 

A = Z^e^Ei/STce^kBTa > 1, (2.3) 

for the formation of a well-defined preferential shell. Here the molecular volumes of the 
two componenta are assumed to be of the same order vo and we define 

e, = de{<P)/d<p, (2.4) 

where the derivative is at fixed T and p. In these relations, we have estimated the 
binding energy eb of a polar molecule to an ion to be Afc^T or Aq^bT. Since eb is a 
microscopic energy, these estimations should be very crude (as well as the Born formula). 
In this theory, the strong solvation regime is realized for ^ 1 or for A ^ 1. If the 
infiuence on the hydrogen bonding network is neglected, the radius -Rgheii ^ solvation 
shell grows as A^^^a or as A^^^a from the decay of the electrostatic energy density oc r~^. 
Here it is worth noting that the shell radius -Rgheii pl^-ys the role of the hydrodynamic 
radius R^^ in the diffusion constant Di = ksT /QirriR^ with rj being the shear viscosity. 
As a result, ions with smaller molecular radii (~ R\on) have larger hydrodynamic radii, 
which is explained from the relation R\ = Rl^^^ [39]. For example, Li"*" ions are smaller 
than Na^ ions, but the diffusion constant of Li"*" is smaller than that of Na"*" in water. 
In electrochemistry, much attention has been paid to the ion distribution and the 



electric potential difference across a liquid-liquid interface 39 , 40 . Remarkably, water 



molecules solvating a hydrophilic ion in the water-rich region remain attached to the 
ion even when the ion crosses an interface and moves into the oil-rich region. Osakai et 
al |4l] measured the amount of water molecules extracted together with hydrophilic ions 
in a nitrobenzene(NB)-rich phase with a small water composition (~ 0.168M). In NB- 
water at room temperatures, the number of coextracted water molecules in a NB-rich 
phase was estimated to be 4 for Na"*", 6 for Li"*", and 15 for Ca^'^. Furthermore, using 
proton NMR spectroscopy, Osakai et al. |42] studied successive formation of complex 
structures of anions (such as Cl~ and Br~) and water molecules by gradually increasing 
the water composition in NB. 

These experiments demonstrte that the solvation shell composed of water molecules 
is stable even when the ambient water composition is much decreased. This is a 
remarkable result suggesting a number of important consequences. For example, with a 
small amount of hydrophilic ions in oil, let us gradually increase the water composition 
from zero. If the binding energy tb of a water molecule to such an ion is much larger 
than A;_bT, we may argue that a hydration shell is formed for [2lj 

> 0soi ~ exp(-eb//i;BT), (2.5) 

where the crossover composition (p^oi is very small for strongly hydrophilic ions. The 
solvation chemical potential iJ^l^ii^) hydrophilic ions in water-oil should decrease 
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considerably in the narrow dilute range < < (j)so\- Its variation should be milder 
in the range (2.5). The solbility of such ions should also increase abruptly in the same 
narrow range of with addition of water to oil. Therefore, solubility measurements of 
hydrophilic ions are informative at small water compositions. 



2.2. Gibbs transfer free energy and solvation chemical potentials 

We consider a liquid-liquid interface between a polar phase a and a less polar phase /3 
with bulk compositions (pa and (^^, across which there arises a difference in /i'^j due to 
its composition dependence: 

= - & (2.6) 
where /i*^ {K = a, (3) are the bulk values of the solvent chemical potential of species 



i in the two phases. In electrochemistry 40,41 , the difference of the solvation free 
energies AG^^ between two phases have been called the standard Gibbs transfer energy. 
(Since AG^^ is usually measured in units of kJ per mole, dividing it by the Avogadro 
number gives A/i^^ per ion). It is well-known that if there are differences among A/i^^ 
(z = 1, 2, • • •), an electric double layer emerges at the interface, giving rise to an electric 
potential jump A$ = $a — ^p-, across the interface in equilibrium (called the Galvani 
potential difference). Similar potential differences also appear at liquid-solid interfaces 
(electrodes) [39] . It is also worth noting that there can be an electric potential difference 
between two liquids due to orientation of molecular dipoles even without added ions, 
for which see Ref.[43] for water-hexane. 

To examine the above effect, let us suppose only two species of ions (z = 1, 2) with 
charges Zie and {Zi > and Z2 < 0). At sufficiently low ion densities, the ion 



chemical potentials are expressed as 39 , 40] 

fii = ksT In(niA^) + Z^e* + fil,{<j)), (2.7) 

for 2 = 1,2. The Aj is the thermal de Broglie length (but is an irrelevant constant in the 
following). From the charge neutrality in the bulk regions, the bulk ion densities nia, 
^1/3, n2a, and satisfy 

ZiUia + Z2n2a = 0, Zinii3 + ^2722/3 = 0. (2.8) 

The continuity of /ij across the interface yields 

kBT ln(n,„/n,^) + Z,eA$ + A/xj,^ = 0, (2.9) 



for i = 1, 2. Then the Galvani potential difference is expressed as 11,40 



A* = ^§1-^, (2,10) 

e(Zi + IZ2I) 

The ion densities in the bulk two phases are related by 

In f ^) = In = + (2.11) 

ynipJ \n2pJ {Zi + \Z2\)kBT ' ^ ' 
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Let us take the z axis perpendicular to a liquid-liquid interface. While the interface 
thickness is given by the correlation length ^, the potential ^{z) changes on the scale 
of the Debye-Hiickel screening length, in phase a and in phase /3. As a result, 
$ changes from to $^ on the spatial scale of + k^^, which becomes very long 
as the ion densities determined by equation (2.11) become very small in one of the two 



phases 11 . Remarkably, A$ in equation (2.10) is independent of the ion densities. It 
is typically of order 10/cBT/e(= 258mV) for A/i^^ 7^ ^f^a/s- It vanishes for symmetric 
ion pairs with A/i^^ = A/i^^. As equation (2.11) indicates, the ion densities in the two 
phases are very different in most cases. For example, if A/i^^/fc^T = A/i^^/fc^T = 10 in 
the monovalent case, the common ratio nip/nia = n2p/n2a becomes e~^° = 4.5 x 10~^. 
If is macroscopic (as in air in contact with salted water), A$ may not be observable. 
In the literature, data of AG^^ on water-nitrobenzene at room temperatures are 



available 40 41 , where water and NB are strongly segregated but a considerable amount 
water is present in the NB-rich phase (0.168M). For aqueous mixtures with a being the 
water-rich phase, AG^^ is negative for hydrophilic ions and positive for hydrophobic 
ions. For water-nitrobenzene |40j, we have A/i^^/fc^T = -13.6 for Na^, -15.3 for Li+, 
-26.9 for Ca^"*", -11.3 for Br~, and -7.46 for I~ as examples of hydrophilic ions, while we 
have A/i^^/fc^T = 14.4 for hydrophobic BPhJ. In these experiments, a hydration shell 
should have been formed around each hydrophilic ion even in the water-poor phase (3. 
Equation (2.5) then suggests (pp > 0soi- Though more data are needed, those of AG^^ 
for NB- water indicates the composition-dependence of filoi{4>) quantitatively, which is 
very strong even after the formation of solvation shells. 

3. Ginzburg-Landau Theory 

We will first give a short summary of our previous results in the Ginzburg-Landau 
scheme in the presence of the preferential solvation in Subsections 3.1-3.3. New results 
will then follow on the structure factors among ions, the surface tension behavior with 
antagonistic salt, and the phase diagram. We treat two species of monovalent ions. 
However, for three species of ions, we obtain more complex ionic distributions near 



liquid-liquid interfaces 12 , 30 and solid surfaces 126 



3.1. Electrostatic and solvation interactions 

We consider a polar binary mixture (water-oil) containing a small amount of monovalent 
salt {Zi = 1 and Z2 = —1). The ions are sufficiently dilute and their volume fraction is 
negligible. Hereafter the Boltzmann constant will be set equal to unity. Neglecting the 



image interaction, we assume the free energy density 10, 11 



f = f,{<P,T) + ^Cm' + ^E' 



+ T 



J2[HniX^)-^-9i<p]nr. (3.1) 
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The first term /o is the chemical part. In our simulation we use the Bragg- WiUiams 
form, 

^/o = 01n0 + (1 - 0) ln(l - 0) + - 0), (3.2) 

where Vq is the solvent molecular volume common to the two components and x = x{T) 
depends on T (at constant pressure). Here x = 2 at the solvent critical temperature 
T = Tc without ions and x — 2 oc T — for small T — Tc. The second term is the 
gradient part with C being a constant. The third term is the electrostatic free energy, 
where E = — V$ is the electric field and the potential $ is produced by the charge 
density p = e{ni — ^2) via the Poisson equation, 

- V ■ e(0)V$ = 47rp. (3.3) 



The dielectric constant e{(j)) can depend on the composition 38 . In our previous work 
the linear composition dependence 

e(0)= £0 + ^10 (3.4) 

has been assumed, where Eq and Si are constants. The last term in the right hand side 
of equation (3.1) consists of the entropic part and the solvation part of the ions, with Aj 
being the thermal de Broglie length. The solvation terms (oc gi) follow if /isoi(0) depend 
on linearly as 

/^soi(</') = const. - T^i0, (3.5) 

where the constant term gives rise to a contribution linear in rij in / and is irrelevant, 
while the second term (oc Qi) yields the solvation coupling in /. In this approximation, 
the parameter Qi represents the strength of the preferential solvation of ion species i. 
The difference of the solvation chemical potential in two-phase coexistence in equation 
(2.6) is given by 

A/i^„i = TgA<P, (3.6) 

where A0 is the composition difference. Thus Qi > for hydrophilic ions and Qi < for 
hydrophobic ions. The discussion in the subsection 2.2 indicates Qi ~ 15 for Na~^ ions 
and Qi ^ —15 for BPh4 in water-nitrobenzene at 300K. 

The linear form equation (3.5) is adopted for the mathematical simplicity and 
should not be taken too seriously. Moreover, around equation (2.5), we predict a steep 
drop of yUsoi(0) in the range < < 0soi for hydrophilic ions with formation of a 
hydration shell. Thus equation (3.5) is applicable for > 0soi. In addition, in our 
previous papers 11 , 12 , the image interaction was also included in the free energy, but 
it is neglected in this paper for simplicity. 

For a monovalent antagonistic salt, figure 2 gives equilibrium interface profiles with 
gi = —g2 = 10 at X = 3 (see Ref. [11] for more details), where we can see a marked 
electric double layer. The bulk cation densities in the two sides, nia and nip, coincide 
from equation (2.8) and is set equal to 2 x 10^'^t'g"^. The interface thickness ^ is of 
order 5a, while ^{z) changes on the scale of ~ ~ 10a. The surface tension 
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Figure 2. Normalized potential $(z)/A$ and composition 4){z) (upper panel), and 
normalized ion densities VQni{z) and wo't-2(^) (lower panel) (taken from Ref.[ll]), 
where X = 3,(7i =52 = 10, and Wo'^iq = t'o^-i/S = 2 x 10"''^. Here A = A means 
jasoT = 8.5. Hydrophilic cations (hydrophobic anions) tend to be on the left (right) 
near the interface, resulting in a large electric double layer. The asymmetry of the two 
curves of ni and 712 is due to the 0-dependence of e{<f>) (ei/eo = 4/3 here). 





change due to salt is Aa = a — = — 0.041Ta~^ (where ctq = 0.497T/a^ without salt) 
and the surface adsorption of ions is F = 0.014a~^ (see equation (3.33)). If gi = —g2, 
the electric double layer is increasingly enhanced with further increasing gi and/or nia, 
eventually leading to vanishing of the surface tension (see figure 3). We can understand 
these features by solving the nonlinear Poisson-Boltzmann equation 12 . See also a 



similar recent calculation of the interface profiles 19 . The ionic interfacial distribution 



was measured for a ternary, antagonistic salt in water-NB in the experiment by Luo et 



al. 30 and was theoretically explained in our previous paper 12 



3.2. Dynamic equations 



We present the dynamic equations for 0, ni, 77-2, and the velocity field v 16,44 . The 
fluid is incompressible and isothermal, so we assume V-i> = and treat the mass density 
Po and the temperature T as constants. Then cf) and rii obey 



dcj) 
di 



+ V- 



T 66' 



(3.7) 
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dt 



+ V 



[niv) = —V ■ riiV — 

I OUi 



Vui =F T^UiE 



giUiVcp 



(3i 



The kinetic coefficient Lq and the ion diffusion constants Di and D2 are constants. The 
symbol ± in the second hne of equation (3.8) denotes — for the cations {i = 1) and + 
for the anions {i = 2). The momentum equation is expressed as 
dv 



Po 



dt 



-Vpi + T]V^v - V ■ JT, 



where pi ensures V ■ v = 0, rj is the shear viscosity, and 77 
the reversible stress tensor of the form, 

n„/3 = TCV„0V/30 - ^E^Ep. 



(3.9) 

{n„/3} {a,/3 = x,y, z) is 



(3.10) 



with = d/dxo,. It follows the relation V ■ J7 = (j)V{6F/A(j)) + Y,.niV{6F/6ni), 
which is needed for the self-consistency of our model dynamic equations [16} [44] 

In numerically integrating the dynamic equations, we may use the following two 
approximations p!6|. First, in the Stokes approximation, we set the right hand side of 
equation (3.9) equal to zero. Then v is expressed in terms of the Oseen tensor inversely 
proportional to the viscosity r] 44 . Second, in the quasi-static approximation, we 



assume that the characteristic ion diffusion time is much faster than the the relaxation 
time of the composition. This is justified near the solvent critical point and in describing 
slow domain growth. In such cases, the ion distributions are expressed in terms of (j) 



and $ in the modified Poisson-Boltzmann relation 11 
n° exp[^i0^e$/T]. 



(3.11) 



This form was also used in analysis of experimental data 30 and in that of molecular 
dynamics simulation 25 . The coefficients n° are determined from the conservation of 
the ion numbers, 



(Ui) 



drriAr) 



no, 



(3.12) 



where (■ ■ ■) = V^'^ J dr{- ■ •) denotes the space average with V being the cell volume. 
The average density = (ni) = (^2) is a given constant. It will be measured in units 
of Vq^ or the normalized density VqUq will be used with vq appearing in equation (3.2). 



3.3. Structure factors and effective ion-ion interaction in one-phase states 

We suppose small deviations 6(j){r) = (f){r) — (0) and 5nj(r) = ni{r) — uq of the 
composition and the ion densities in one-phase states. Their equilibrium distributions 
are treated to be Gaussian in the mean field theory. Then, from the free energy (equation 
(3.1)), we may readily calculate the following structure factors, 

C{q) = MjeW (3.13) 
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where (f)q, n^g, and pg are the Fourier components of the composition (j), the ion densities 
Hi, and the charge density p = e{ni — ■^2) with wave vector q and (■ ■ ■)e denotes the 
thermal average. It follows the relation, 

C{q) = Gu{q) + G22{q) ~ 20^- (3-14) 

The 0-dependence of the dielectric constant e is negligible for small fluctuations. We 
introduce the Bjerrum length ^B and the Debye wave number /t by 

£g = eVeT, K = (Svr^Br^o)^^^- (3.15) 

After elimination the ion fluctuations, the second order contribution of the free 
energy is written as = T 1 0^^/21^ S'(g). The inverse of S{q) is written as 

,2,^2 



10 



11 



^ = r - (^71 + g,?"^ + Cq' \l - , (3.16) 

where f is related to the second derivative of /o and use of equation (3.2) yields 

f = f^/T = v^'[l/<^{l - 0) - 2x]. (3.17) 

The parameter 7p represents asymmetry of the solvation of the two ion species and is 
defined by 

7p = (167rC£B)-i/2|(7i-(?2| (3.18) 

Remarkably, 7p is independent of the salt density. Note that the structure factor 
for weakly ionized polyelectrolytes has the same form [2}[3], where the solvation has 
been neglected, however. Recently, we have presented a generalized form of the 
structure factor for polyelectrolytes including the solvation interaction and the ionization 
fluctuations [20] . 

We explain implications of the structure factor (equation (3.16)). (i) The second 
term gives rise to a shift of the spinodal curve. For example, if the cations and anions 
are hydrophilic and f/i ~ 5^2 ~ 15, the shift term is of order — 500no and its magnitude 
can be appreciable even for vqUq <^ 1. It is well-known that the coexistence curve and 
the spinodal curve of aqueous mixtures are much shifted even by a small amount of 
hydrophilic ions such as Na"*" and Br~ |45j. (ii) On the other hand, 7p can be increased 
for antagonistic salt |To|[Tl||3l]. From the last term in equation (3.16) a Lifshitz point 
appears at 7p = 1. (In figure 2, 7p = 5/4v^ < 1). That is, while S{q) is maximum at 
g = for 7p < 1, it exhibits a peak at an intermediate wave number q^ = (7p — 1)^/^k 
for 7p > 1. The peak height is given by S'(gm) = l/(^ ~ ''^m), where 

rm = (91 + g2?n^/2 - C(7p - 1)'^^. (3.19) 

For f < Tm, mesophase formation should take place (see figure 4). For f > rm, two- 
phase states with a planar interface can be stable or metastable even for 7p > 1. (iii) 
If 7p < 1, the criticality and the macroscopic phase separation can remain in the mean 
field theory. That is, for g ^ k, we obtain the Ornstein-Zernike form, 

S{q) - l/[f - {g, + g^fn,/2 + C(l - ^l)q\ (3.20) 
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The coefficient in front of is reduced by the factor 1 — 7p , which leads to enhancement 
of S{q) in the range < q < k. When k > the correlation length ^ should be 
redefined by 

e = C{l-^l)/[f-^-{g, + g,rno]. (3.21) 

Second, retaining the ion densities, we ehminate the composition fluctuations in F 
by setting 

if + Cq^)^g = J29^n^g, (3.22) 

i 

where the terms nonlinear in 6(f) are neglected. We then obtain effective interactions 
among the ions mediated by the composition fluctuations. The resultant free energy of 
ions is written as 



dr Trii ln(ni A ■ ) 

i 

+ ^ ^dr^dr'Y,y^3{\r-r'\)5n,{r)5n,{r'), (3.23) 

where the terms linear in nj are not written explicitly. Here the deviations Srii = ni — {ni) 
need not be very small and the logarithmic terms are not expanded with respect to them. 
The effective interaction potentials Vij{r) are expressed as 



T/.,(r) = ^-^le-/C, (3.24) 
er AnC r 

where Zic and are the ion charges (being ±e in the monovalent case) and 
^ = [C /fY^"^ is the correlation length without ions (different from ^ in equation (3.21)). 

The solvation-induced interaction is effective in the range a ^ r < ^ and can be 
increasingly important on approaching the solvent criticality (for ^ a), (i) It is 
attractive among the ions of the same species [i = j) and dominates over the Coulomb 
repulsion for 

> ATvCiB- (3.25) 

Under this condition there should be a tendency of ion aggregation of the same species. 
(ii)In the antagonistic case {gig2 < 0), the cations and anions can repel one another in 
the range a ^ r < ,^ for 

\gig2\ > AnCiB, (3.26) 

which triggers charge density waves near the solvent criticality. 

The ion structure factors Gij{q) in equation (3.13) can be calculated from Fion in 
equation (3.23). The inverse of the 2x2 matrix Gij{q) is written as 

G'^{q) = 6i,+noV,,{q)/T, (3.27) 
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Figure 3. Surface quantities a, 2ae, '2ag, and (Tq ^ TT + erg in units of Ta~^ vs voUa 
for a monovalent, antagonistic salt with gi = —52 = 10 at x = 2.1. Here 2(jg and 
CTo — rr + Ce nearly coincide to support equation (3.36). 



where Vij{q) is the Fourier transformation of Vij{r) in equation (3.24). In the monovalent 
case, some calculations give desired results, 

Giiiq) ^ , 2 , 2 (91 -92? ^gfu 



Go{q) 



1 + noS{q) 



9i+92- 



2(u + l) 2m + 1 



^o(.) - ^ - (3.29) 



G'i2(g) = 2 + f %)(^7i + ^?2)^-^C7(g), 

^^(g) = ^+-o%)^f^^«^ (3.28) 
where m = / n"^ and 

M + 1/2 _ g2 + ft;2/2 
M + 1 g2 + 

is the structure factor for the cations (or for the anions) in the absence of solvation. The 
solvation parts in equation (3.28) are all proportional to noS'(g), where S{q) is given by 
equation (3.16), so they diverge on approaching the spinodal discussed below equation 
(3.18). The Coulomb interaction suppresses large-scale charge-density fluctuations, so 
C{q) tends to zero as g — )■ 0. 

3.4- Surface tension 

Addition of an antagonistic salt can dramatically reduce the surface tension a of liquid- 



liquid interfaces even at small content 12 . This trend was in fact observed 29 30 . 



Here the behavior of a will be examined in some detail as a continuation of Ref.[ll], 
but the image interaction will be neglected. 

We suppose a planar liquid-liquid interface at z = taking the z axis in its 
perpendicular direction. We impose the boundary condition — 0^ in phase a and 
— )■ 0^ in phase /3 far from the interface and minimize the grand potential ^ = j dru. 
The grand potential density u reads 

f - hep -"^fiiUi, (3.30) 



00 



Solvation Effects in Phase Transitions in Soft Matter 



14 



where / is the free energy density in equation (3.1) and the chemical potentials 
h = 6F/6(j) and yUj = 6F/6ni are homogeneous constants. The relation du/dz = 
2C(f)'(j)" — d{p^)/dz is then derived and is integrated to give 

co = C(p'^ - p^ + cooo. (3.31) 
where 0' = d(f)/dz and 0" = d'^ip/dz'^. Since 0' and p tend to zero far from the interface, 
w(z) tends to a common value as 2; — )■ ±00. The surface tension is written as 

a = j dz[u{z) - Woo] = j dz[C(t)'^ - (3.32) 

In our previous work (ll ]l2 , we expanded cr with respect to the salt density as 

a^aQ-TT - / dz—E^, (3.33) 

J OTT 

where cxo is the surface tension without solute and T is the surface adsorption of ions. 
In terms of the total ion density n = ni + n2, it is defined by 



f An 

r = y dz[n -Ua - -^{4> - 4>a)], 



(3.34) 



where = niK + n2K {K = a 01 (5) and An = Ua — np, so the integrand tends to zero as 
z — )■ ±00. The expansion (equation (3.33)) is valid up to linear order in the salt density. 



For dilute neutral solute, it becomes the well-known Gibbs formula a = (Tq — TV 46 



The electrostatic contribution to the surface tension has been neglected in the Gibbs 



formula, but it can be crucial for antagonistic salt 12 and for ionic surfactant 22 



We have numerically examined these relations. To this end, we introduce the areal 
densities of the electrostatic energy and the gradient free energy as 

cTe = / dzeEySTc, (Tg= I dzC(j)'^/2. (3.35) 



Then equations (3.32) and (3.33) yield a = 2ag — 2ae and 

2ag ^ CTo - + de. (3.36) 

In figure 3, we show the surface quantities 2ae, 2ag, a = 2ag — 2ae, and the combination 
(To — TT + (Te as functions of the bulk ion density n^, where gi = —g2 = 10, x = 2.1, 
and 7p = 3.455. Here = Uf^ for gi = —g2- The surface tension vanishes at 
vo^a = 2.5 X 10^^, where ag = Oe- The approximate relation (equation (3.36)) excellently 
holds here, where its deviation is at most of order 10~^T/a^. 



3.5. Phase diagram with antagonistic salt 

With addition of antagonistic salt, the consolute criticality of binary mixtures disappears 
due to emergence of mesophases. However, we have not yet fully understood the phase 
behavior even in the mean field theory. So far, we have found (i) instability of one- 
phase states for 7p > 1 in the range f < r^, as discussed around equation (3.19), and 
(ii) vanishing of the surface tension a with increasing the salt density as in figure 3 or 
on approaching the solvent criticality. Thus, in figure 4, mesophases should be formed 
in the area with x < 2.14 and f < in the (p-x plane. Here gi = —g2 = 15 at 
VqUi = 4 X 10""^. 
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Figure 4. Phase diagram for gi = —g^ = 15, VqUi = 4 x 10"'', and 7p = 1.727 
in the x-4> plane. Mesophases are reahzed in the region (in green) where x < 2.14 
and f > Tm- Macroscopic phase separation occurs in the two-phase region above the 
mesophase region. 




Figure 5. Structure factor S{q, t) vs q at t = lO^io of a binary mixture at {(p) — 0.4 
and novo = 0.002 for (a)x = 2.0, (b)2.1, (c)2.2, (d)2.3, and (e)2.4. 



4. Simulations in Two Dimensions with Antagonistic Salt 

We already presented 2D simulation results of phase separation of a binary mixture 
with a monovalent, antagonistic salt at the critical composition (0) = 0.5 [Tg]. Here 
we present further results in two dimensions at (0) = 0.5 and 0.4. The composition 
patterns are analogous to those well-known for block copolymers, surfactant systems, 



and magnetic fluids 44,47,48 . In our ionic case, however, we shall encounter some 
unique features in the droplet patterns. We stress that three dimensional simulations 
are needed to really understand the phase ordering in this system. For example, we 



should reproduce the observed onion structure 32 . 
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Figure 6. Characteristic domain size 27r/qm(t) vs i of a binary mixture at (0) = 0.5 
(top) and 0.4 (bottom) for n^VQ = 0.002 for (a)x = 2.0, (b)2.1, (c)2.2, (d)2.3, and 
(e)2.4. For {(f) — 0.5, the domain structure is nearly pinned at fong times for these x- 
For {(()) = 0.4, the coarsening is nearly stopped for x = 2.0 and 2.1, while it is much 
slowed down for the larger x- As a guide, a line with slope 1/3 is written (dotted line). 



(a) 1 = 1.9 {b) X=2. l {c)y=2.2 (d) 7=2.4 




Figure 7. Composition patterns of a binary mixture at (0) ~ 0.5 containing an 
antagonistic salt with vqUq = 0.02 (top) and 0.03 (bottom) &t t — lO'^io, where 
(a)x — 1.9, (b)2.1, (c)2.2, and (d) 2.4. These patterns are frozen in time. 
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Figure 8. Composition patterns of a binary mixture at (0) = 0.4 containing an 
antagonistic salt with vqUq — 0.02 (top) and 0.03 (bottom) at t = lO^'io, where 
(a)x = 2.0, (b)2.1, (c)2.2, and (d) 2.3. Patterns are fronzen for (a) and (b) and 
are very slowly evolving for (c) and (d). Size dispersity of droplets is conspicuous. 



Background of simulation 

We vary the interaction parameter x ^i-nd set the average ion density no = {rti) = (77,2) 
equal to 0.002t;(^^ or 0.003f(j'^. The other parameters are the same as in our previous 
work [16j. Namely, we set C = a~^, gi = —g2 = 15, ei = 0, and Eb = 3a. Then we 
obtain 7p = 2.44 from equation (3.18). These values of gi and g2 are realistic in view 
of the data of the Gibbs transfer free energy, as discussed below equation (2.11). On 
a 256 X 256 square lattice, we integrated the dynamic equation (3.7) using the Bragg- 
Williams free energy density fo in equation (3.1), the Stokes approximation for the 
velocity field v, and the modified Poisson-Boltzmann expressions (equation (3.12)) for 
the ion densities rii and n2. At t = 0, we started with the initial condition (f){r, 0) = (0) 

with small random numbers superimposed. Subsequently, phase ordering took place for 

1/3 

the parameter values given below. The space mesh size is a = Vq and the time step of 
integration is O.Olto, where to is related to the kinetic coefficient Lq in equation (3.7). 
and the viscosity 77 in equation (3.9) as 

to = aV^o, Lo = 0.16a^T/ri. (4.1) 

Time-evolution of the composition is meassured in units of to- In late stages it is much 
slower than to- 

4-2. Slow coarsening and patterns 

In figure 5, we plot the (angle-averaged) structure factor S{q,t) of the composition vs 
q for {(f)) = 0.4 at t = lO'^to- It exhibits a sharp peak at g ~ qm.{t)- In our previous 
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<(j)>= 0.4 




Figure 9. Composition cation density ni, anion density 712, and potential $ (from 
top to bottom) for x = 2-1 (left) and x = 2.4 (right) in a binary mixture at {(/)) = 0.4 
and uoVq — 0.002. Snapshots in a 1/4 of the total system are displayed. 



work 16 , a similar plot of S{q,t) was given for (0) = 0.5. For each simulation run, we 
calculate the characteristic wave number qm{t) by 

qm{t) = J2\qmVj2\^,\\ (4.2) 

where (j)q is the Fourier component of (l){r,t). In figure 6, we plot the characteristic 
domain size 2Tr/qm{t) vs t. For all x (0) = 0.5 and for x =2.0 and 2.1 at (0) = 0.4, 
the coarsening is nearly stopped at long times. For x ^ 2.2 at (0) = 0.4, it is much 
slowed down at long times and, if we write ~ t"^ , the exponent 7 increases from 
zero up to about 0.1 at x = 2.4. Note that the typical domain size grows as t^^^ for 
conserved order parameter without the hydrodynamic interaction and is faster with the 



hydrodynamic interaction 44 ,49 



We next show composition patterns at t = lO^to- From figure 6, they are frozen or 
very slowly evolving. In figure 7, they are stripe-like and bicontinuous for (0) = 0.5. In 
figure 8, they are droplet-like for (0) = 0.4. The gradation of each snapshot represents 
the local value of 0(r,t). For (0) = 0.5 (or 0.4), the variance [((0 — (0))^)]"'^''^ is about 
0.16 (or 0.13) at X = 2 and increases to 0.33 (or 0.31) at x = 2.4. Note that no phase 
separation occurs for x ^ 2 without ions. We recognize that the patterns are very 
sensitive to the values of x and fo^o- A marked feature is that the largest domain 
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Figure 10. Effective surface tension in equation (4.3) vs t for {(f)) = 0.5 (top) 
and 0.4 (bottom). It tends to zero for (0) — 0.5 but remains positive for {(f) = 0.4. 



size increases with increasing x ^i^nd decreases with increasing uq. More details are 
as follows, (i) At (0) = 0.5 in figure 7, the patterns are lamellar-like for x = 1-9) 
while they resemble those in systems with competing interactions 47,48]. (ii) In the 
droplet patterns at {(f)) = 0.4, water-rich droplets in a percolated oil-rich region have 
considerable size dispersity. This is possible only when even small droplets have long 
life times. (Without ions, small droplets quickly disappear due to the surface tension 
effect.) Large and small droplets segregate exhibiting polycrystal-like medium-range 
order for x = 2, while larger droplets are growing and smaller ones are diminishing very 
slowly in random configurations for x ^ 2.2. 

In figure 9, at {(f)) = 0.4, we display the composition (f){r, t), the ion densities ni{r, t) 
and n2(r,t), and the potential $(r,t) at t = lO'^to for the two cases x = 2.1 and 2.4. 



Similar snapshots at (0) = 0.5 can be found in our previous paper 16 . The typical 



droplet radii are a few times larger for x = 2.4 than for x = 2.1, but small droplets 
(with diameter~ 5a) still remain in the two cases giving rise to sharp peaks in the cation 
distribution. Due to the strong selective solvation, the cations are within the droplets 
and the anions are outside them. For x = 2.1, the typical droplet size is small and the 
cation density is highest at the droplet centers. For x = 2.4, the cations and the anions 
are mostly localized in the electric double layers at the interfaces. In these cases, the 
water-rich droplets are positively charged and are repelling one another. However, we 
have not yet examined how this repulsion can infiuence the coarsening process. 
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4-3. Effective surface tension 

In figure 10, we show the time-development of the effective surface tension defined by 

aeff = ^ 1 driCm' - (4.3) 

where S is the total interface length calculated separately. This quantity becomes the 
usual surface tension for a planar interface in equilibrium from equation (3.32). (In 



our previous work 16 , the space averages of the gradient free energy density and the 
electrostatic energy density were separately plotted for (0) = 0.5 in figure 10 there.) 
We recognize that a^s tends to zero for bicontinuous domain morphology but remains 
positive for droplet morphology. It is rather surprising that the coarsening in the droplet 
case is much slowed down even for positive aes- To explain these results theoretically, 
we need more systematic analysis of the problem. 



5. Summary and remarks 

In this work, we have demonstrated the crucial role of the preferential solvation of ions 
in phase transitions of polar binary mixtures. Together with a short summary of our 
previous work, we have presented some new results, (i) In equation (3.28), we have 
calculated the structure factors among the ion density fluctuations, (ii) In figure 3, 
we have examined how the surface tension tends to zero with increasing the density of 
antagonistic salt, (iii) In figure 4, we have shown a phase diagram of a binary mixtures 
containing antagonistic salt, where mesophases are realized near the solvent criticality. 
Note that the original criticality without ions disappears when antagonistic salt with 
7p > 1 is added, (vi) In figures 5-10, we have numerically examined the phase ordering 
of a binary mixture containing antagonistic salt for four values of x in the bicontinuous 
case at (0) = 0.5 and in the droplet case at (0) = 0.4. The coarsening is markedly 
slowed down with decreasing x to 2 and/or with increasing the ion density no- The 
size dispersity of water-rich droplets is significant, because even small droplets have 
very long life times. In figure 9, we have found that water- rich droplets are positively 
charged with the anions in the percolated oil- rich region. In figure 10, the effective 
surface tension a^s vs t is plotted for (0) = 0.5 and 0.4. It tends to zero for (0) = 0.5 
but remains positive for (0) = 0.4. 

In figures 7 and 8, the domain configurations exhibit some mesoscopic order 
(lamellar-like or polycrystal-like) for x close to 2, while they remain disordered for 
larger x in our simulation time. In this paper, we cannot conclude whether or not the 
domain growth has ultimately stopped or not. If the thermal noise is present, the system 
might be slowly approaching lamellar or crystalline states at very long times. 

In our theory, the molecular volumes of the two components are commonly vq- 
However, the molecular volumes of D2O and 3MP (the inverse densities of the pure 
components) are 28 and 168 A^, respectively, for example. Moreover, the coefficient 
C of the gradient free energy is an arbitrary constant. Therefore, our theory remains 
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very qualitative. It is not clear how C can be determined in aqueous mixtures where 
the hydrogen bonding is of primary importance. For pure water, the observed surface 
tension outside the critical region 1 — T/Tc ^0.1 fairly agrees with the calculated 
surface tension from the van der Waals model with the gradient free energy density 



10Ta^|Vnp/2 52,53 , where n is the density and a = 3.lA is the van der Waals radius. 



There can be a number of dynamical effects in the presence of antagonistic 



salt, including multi-exponential decays of the dynamic light scattering intensity 16 



and singular rheological behavior with appearance of mesophases |44;|. We propose 
experiments of the surface instability taking place for negative surface tension. That 
is, a water-oil interface becomes unstable with addition of antagonistic salt above a 
critical amount (as suggested by figure 3) or on approaching the solvent criticality with 
a fixed density of antagonistic salt. In the introduction, we have already mentioned 
the experiments of spontaneous emulsification at a water-oil interface with antagonistic 



salt 33 34 



Though this paper has mostly treated antagonistic salt, intriguing solvation effects 
can well be predicted for aqueous mixtures containing usual hydrophilic ion pairs. 
For example, many groups have detected long-lived heterogeneities by dynamic light 
scattering in one-phase states of various aqueous mixtures containing hydrophilic salt 



50 . Recently, we have shown that preferential solvation can stabilize water-rich 



domains enriched with ions even for x < 2 when water is the minority component 51 



We also mention the problem of preferential wetting of mixture solvents on ionizable 
surfaces, which can be either planar, spherical, or cylindrical. Around an ionizable rod, 
we already predicted a first-order wetting transition with discontinuous jumps in the 



preferential adsorption and the degree of ionization 21 . In such situations, we should 



take into account the solvation interaction among the charges (consisting of the ionized 
units on the surface and the counterions) and the polar molecules. 

Finally, it is worth noting that ion-induced nucleation has long been investigated 



in the literature 54 , 55 . It is well-known that water droplets can easily be produced 
around hydrophilic ions in metastable gas mixtures containing water vapor. Here ions 
play the role of nucleation seeds on which hydration-induced condensation occurs. A 



Ginzburg-Landau approach to this problem in the line of Ref. [9] was also presented 53 
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